Unified optical-model approach to low-energy antiproton 
annihilation on nuclei and to antiprotonic atoms 



C.J. Batty", E. Friedman^ A. Gal'' 
"Rutherford Appleton Laboratory, Chilton, Didcot, Oxon 0X11 OQX, UK 
^Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 



Abstract 

A successful unified description of p nuclear interactions near ii^ = is 
achieved using a p optical potential within a folding model, V^pt ^ v * p, 
where a pp potential v is folded with the nuclear density p. The potential 
V fits very well the measured pp annihilation cross sections at low energies 
{pl < 200 MeV/c) and the Is and 2p spin- averaged level shifts and widths 
for the pH atom. The density-folded optical potential V^pt reproduces satis- 
factorily the strong-interaction level shifts and widths over the entire periodic 
table, for A > 10, as well as the few low energy p annihilation cross sections 
measured on Ne. Both v and Vopt are found to be highly absorptive, which 
leads to a saturation of reaction cross sections in hydrogen and on nuclei. 
Predictions are made for p annihilation cross sections over the entire peri- 
odic table at these very low energies and the systematics of the calculated 
cross sections as function of A, Z and E is discussed and explained in terms 
of a Coulomb-modified strong-absorption model. Finally, optical potentials 
which fit simultaneously low-energy p— ^He observables for E < as well as 
for E > are used to assess the reliability of extracting Coulomb modified p 
nuclear scattering lengths directly from the data. The relationship between 
different kinds of scattering lengths is discussed and previously published sys- 
tematics of the p nuclear scattering lengths is updated. 
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I. INTRODUCTION 



The present paper is motivated by the recent pubhcations of experimental annihilation 
cross sections for antiprotons on several targets at very low energies [pi < 100 MeV/c) 
Some unexpected features observed in the dependence of these cross sections on the 
mass (A) and charge (Z) numbers have been very recently explained by us 0] as arising 
from saturation of the p annihilation cross sections due to the strong absorption that is 
characteristic of the p-nucleus interaction. Using an optical model approach we have shown 
that the E = borderline between p atoms and the scattering regime can be crossed using 
the same p-nucleus potential, i.e. one can get good fits to experimental results both for 
p atoms and for low energy annihilation cross sections with the same potential, at least 
within a given mass region of the periodic table. Such potentials are strongly absorptive, 
which leads to a remarkable saturation of the total reaction cross section with increasing 
A. Strong absorption has very recently been shown P,^] to lead also to saturation of the 
widths of p atomic states and to the prediction of relatively narrow deeply bound p atomic 
states. The close analogy between bound-state widths and total reaction cross sections is 
best demonstrated by observing the corresponding expressions, assuming for simplicity a 
Schrodinger-type equation. The width F of an atomic level is then given in terms of the 
optical potential Vopt by: 

r ^ /ImKpt(r)|^(r)prfr 

2 /|V(r)|2dr ■ ^ ^ 

Here V^(r) is the p full atomic wavefunction. The corresponding expression for the total 
reaction cross section at positive energies is 

^/? = -^/lmKpt(r)|V'(r)prfr . (2) 

Here ipi'"^) is the p-nucleus elastic scattering wavefunction and v is the cm. velocity. We 
note that these expressions involve no approximation. 

In the present work we extend the earlier work by presenting a full analysis of the 
dependence of calculated annihilation cross sections on A, on Z and on the energy up 
to 16 MeV {pl = 175 MeV). Beginning with the pp system, we show that the available 
experimental results, including the strong-interaction spin-averaged shift and widths of the 
pH atom, are very well accounted for by a potential approach. Moving over to p-nuclear 
systems, we observe that the strong-interaction level shifts and widths in p atoms are well 
reproduced for A > 10, over the entire periodic table, by folding the above pp potential v 
with the nuclear density and renormalizing slightly its strength, to obtain a p-nucleus optical 
potential. Using this density-folded optical potential Vopt, we study the systematics of the 
predicted cross sections over the entire periodic table and show it to result naturally from 
the absorptive properties of the interaction under conditions of Coulomb focussing. The 
optical potentials are then used to derive p-nucleus s-wave scattering lengths. We study the 
A dependence of these scattering lengths and compare them with those extracted directly 
from the low-energy data using approximation methods 

The paper is organized as follows. In Section |I| we define the optical potential and 
the various scattering amplitudes which are used throughout the present work. Section 
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ri| deals with the pp system and in Section |^ we discuss the p interaction with nuclei. 
Section ^ is devoted to the p - ^He system which is used for closely examining the scattering 
length approximation and Section [V^ extends the discussion to heavier systems. The work 
is summarized and concluded in Section IVT] . 



II. OPTICAL POTENTIALS AND SCATTERING AMPLITUDES 

The interaction of low energy antiprotons with nuclei, as well as the interaction of an- 
tiprotons bound in an atomic system, is described in this work by the conventional 'tp' 
potential [|| 

2/iKpt(r) = -47r(l + ^^^)VW , (3) 

where m is the nucleon mass, /i is the p-nucleus reduced mass, 60 is an 'effective scattering 
length' complex parameter obtained from fits to the data and p{r) is the nuclear density 
distribution normalized to A. The density p(r) may also include the effect of folding in a 
finite-range two-body interaction. The factor {A — 1)/A in Eq. (|^) is usually omitted in 
large A atomic studies as we do also here (Sect. |V|), but is retained of course in the 
discussion of the pp system in Sect. |IT1| and the p-He system in Sect. 0. 

It is customary to describe the interaction of bound hadrons with the nucleus by the 
Klein-Gordon (KG) equation of the form: 



V2-2/i(5 + Kpt + K) + (K + 5)2]v; = (n = c=l) , (4) 

where B is the complex binding energy and Vc is the finite-size Coulomb interaction of the 
hadron with the nucleus, including vacuum-polarization terms. Equation (^ assumes that 
Kpt behaves as a Lorentz scalar. The use of the KG equation for p atoms is justified as long 
as spin effects are negligible and one is interested in (2j + 1) averaging. As we are interested 
in this work in crossing the E=Q borderline between the bound atomic system and the low 
energy scattering regime, we use the same wave equation also for positive energies, namely 

V^ + e- (242k - K') - 2/iKpt] ^ = Q {h = c=l) (5) 

where k and are the wave number and reduced energy respectively in the cm. system, 
and (^l-^)"^ = + E^^ in terms of the cm. total energies for the projectile and target, 
respectively. Equation (^ is the Lorentz scalar version of the KG equation used in Ref. 
for studying interaction with nuclei. Near E = {] both equations yield practically the 
same numerical results. Furthermore, we have verified that using the Schrodinger equation 
instead of the KG equation leads to calculated reaction cross sections that differ by 0.5% or 
less. 



The scattering amplitude due to the wave equation (|g) can be written in the form |]10 

/(^) = /W(^) + e2-/(-)(^) , (6) 
where f'^^\6) is the point Coulomb scattering amplitude. 
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and with non-relativistic Coulomb phase shifts given by 

ai = argr(l + 1 - ir]) , f] = , (8) 

with ttB = /{Zjie^) the Bohr radius. If f{6) is defined to have the phase shifts ai + 
with Coulomb-modified phase shifts 6\^^^ due to the short-ranged part of the interaction, 
then the partial- wave expansion of the Coulomb-modified scattering amplitude f'^^'^\0) is 
given by 

/(-)(^) = i^(2i+l)c/^'=)p,(cos^) , ^ e^'t\in5t) , (9) 

subject to the following partial-wave representation of the operator product e^'" f'^^''''{0): 

^2iajisc)^g^ = i ^(2/ + l)e2^'^'Cf '^Pi(cos^) . (10) 

In practice, the Coulomb-modified partial wave amplitudes C;^*'^^ are determined by solv- 
ing the wave equation for each partial wave with point Coulomb boundary conditions at 
sufficiently large distance. The total reaction cross section is then given by 

In addition to the p reaction cross section, which near = is essentially exhausted by 
the total annihilation cross section, we also discuss in this paper the low energy effective- 
range expansion which requires the knowledge of the Coulomb- modified phase shifts Si^'^\ 
These are given in terms of the partial wave amplitudes as follows: 

Im5;^^) = -^ln[(l-2ImC/^'=))2 + 4(ReC/^'=V] , (12) 

[1 - 2ImQ'''^ + 2iReC/''^]e2M - 1 

where the latter expression is particularly useful in order to avoid a possible ambiguity in 
5i^'^^ if it is determined from an alternative expression involving sin(2Re Si^'^^). We note 
that for absorptive potentials T4pt, the Coulomb- modified phase shifts are complex, with 
Im 6^'^'^^ > 0. An alternative expression for in terms of transmission coefficients T; is then 
given by 

^«=^E(2^ + 1)T] , Ti ^ 1 - e-^'-^^'r' , (14) 
which explicitly shows that an—O for real potentials. 
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Finally, the low-energy effective-range expansion for / = 0, in the presence of an attractive 
Coulomb potential, is given by |11,12]: 



where 



Cl iv) k cot 6t^ - A/, (^) = _ 1 + 1 ^ o (^4) (15) 



Co iv) = . _ _2^^ , h{ri) = Re^ {ir]) - In r/ , (16) 



and ^p is the digamma function |T^. For ri>l, or equivalently kas'^ 1, the following limits 
hold: 

Co^(,,l).2vr,, + + (17) 

so that Eq. ( pTSD for the effective range expansion reduces to 

- cot - + ^ (l^B + rt^) k' + O (k^) . (18) 

In the present work we solve numerically the scattering wave equation down to El = 100 keV, 
where the linear dependence of the l.h.s. on k^ is explicitly verified, in order to extract 
reliably the Coulomb modified low energy parameters. We recall, that in the absence of a 
Coulomb potential, the effective range expansion for a short-ranged potential (such as Kpt) 
assumes the form: 

k cot 6o^-— + Kok^ + O (k^) , (19) 
ao 2 ^ ^ 

which can be derived from Eq. ([I5|) upon taking the limit 1] —>■ 0. The Coulomb-modified 

(sc) 

scattering length ag is related to the purely short-ranged scattering length ao by the 
following approximate expression [0,^: 

^ ■ ^ ^111^ + 1-27) , (20) 



a^^'^) ao aB V 2R 

which assumes that R <ti as, where R is the range of the short-ranged potential and 7 
0.5772 ... is the Euler constant. 



III. THE ANTIPROTON-HYDROGEN SYSTEM 



As a first step towards studying the p-nucleus system we tried to fit the data for the 
p hydrogen system with as simple an 'optical potential' v as possible, choosing a Gaussian 
shape common to both real and imaginary parts. In the present application to the pp system 
the 'density' p(r) in Eq. (^) was chosen as proportional to a Gaussian exp(— r^/a|.) and 
normalized to a volume integral of 1. Total reaction cross sections were calculated at 6 
momenta between 38 and 175 MeV/c |1T1,[T6| and the complex parameter Bq was varied in 
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order to fit these recent low-energy data for a given value of ac- This was repeated for 
several values of ac between 1 and 2 fm. In all cases excellent fits to the annihilation cross 
sections were achieved, and consequently the range parameter ac could not be determined 
uniquely from the annihilation cross sections. As we aim at using the same potential on 
both sides of the E = borderline, we turned next to the pB. atom. Using the same type 



of potential as above, excellent agreement between calculation and experiment HT^Ill] was 
achieved by varying bo for any value of ac between 1 and 2 fm. However, requiring that both 
atomic and scattering data be described by the same potential v, we find 0^=1.5 fm and 
bo = —0.15 + zl.8 fm, leading to a combined of 3.8 for the 6 data points for annihilation 
cross sections plus the 3 data points for the atom. The uncertainties are ±0.15, ±0.15 
and ±0.06 fm for ac, Re bo and Im bo, respectively. This potential will be referred to as 
potential (G). A more elaborate potential was previously shown |l9| , p0| to be able to fit the 
low-energy pp annihilation data. In that potential the real part is described by a Saxon- 
Woods (SW) shape with a radius of 1.89 fm and diffuseness parameter of 0.2 fm whereas 
the imaginary part has the SW shape with a radius of only 0.41 fm and the same diffuseness 
of 0.2 fm. Adopting these geometries we repeated the procedure described above for fitting 
the annihilation and atomic data and obtained very good fits with a of 5.8 for the 9 
data points. The parameter bo is 2.85 ± zl6.5 fm, with about ±10% uncertainties, leading 
to potential depths of 46.5 and 7550 MeV for the real and imaginary parts, respectively, in 
very good agreement with Refs. [|^,|^. This potential will be referred to as potential (SW). 
We note that for both potentials (G) and (SW), the imaginary part is considerably stronger 
than the real part, which signals a very strong absorptivity in the pp system. The real 
part of the potential, under such circumstances, plays only a minor role [Q. The imaginary 
parts of the potentials (G) and (SW) have markedly different ranges. The very short range 
of (SW) is qualitatively in agreement with other widely used phenomenological potentials 
pT| , |22[| , whereas the range of (G) is about twice as large as for these potentials. 

Figure |I] shows results for the pp system. The upper part compares calculated reaction 
cross sections for the two potentials with the experimental values, and the excellent agree- 
ment is evident. The lower part shows the ratios of real to imaginary parts of the forward 
Coulomb- modified scattering amplitude f^^^^ [see Eqs. (|]) - (|1^)] for the two potentials, 
which is a quantity of considerable interest: this, so called p parameter, was determined 
at higher p incident momenta to be positive, tending to zero between 200 - 300 MeV/c 



23| . Our calculations give negative values for p at this lower-momentum range, indicating 



a smooth transition to the corresponding ratios for the Coulomb-modified scattering length 
also shown in the lower part, at zero energy, for potentials (G) and (SW). Evidently the 
two potentials, in spite of their different geometry, are very similar in predicting observable 
quantities. 



IV. THE ANTIPROTON-NUCLEUS SYSTEM 



With an established pp interaction potential v, the simplest way of incorporating this 
potential into heavier p-nuclear systems would be to fold it with the nuclear density to 
obtain a p-nucleus potential: Kpt = v * p. Setting aside possible renormalization effects in 
nuclei, this folding procedure also tacitly assumes that at low energy the spin-averaged pn 
interaction is approximately equal to the spin-averaged pp interaction, as supported by the 
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near equality of the imaginary parts of the corresponding scattering lengths and also by 



several phenomenological NN potentials [^,^|- However, recent np low-energy annihilation 



data from the OBELIX collaboration ||2^ suggest that the np interaction, which by charge 
symmetry is equal to the pn interaction, is considerably weaker than the pp interaction. 
Therefore, we keep the strength of the density-folded p-nucleus potential as a free parameter 
in our search, as is discussed below. 

In the energy range covered by the present work the majority of data come from p atoms 
and therefore this naive approach is best tested by applying the resulting potentials to p 
atoms. We emphasize that the p - nucleus potentials obtained from p atoms are not unique. 
For 'macroscopic' nuclear densities of the type discussed in Ref. |^, and for A > 10, a 
Gaussian with a range parameter of about 1.5 fm leads to the best fit. Folding in the pp (G) 
potential with 'macroscopic' nuclear densities, we find that a density-folded optical potential 
of the form @, with bo = —0.1 + il.2 fm, fits all p atoms heavier than boron {A > 10) 
with a of 2.4 per point, whilst the corresponding value for a best-fit phenomenological 
potential of the form Eq. (|) using the same unfolded nuclear densities is 2.7 p. We denote 
this p- nucleus density- folded optical potential by (F). The value of its strength parameter 
6o is 2/3 of the strength bo of the pp potential (G). Note that this is not a best-fit potential 
but merely a renormalized G * p potential, where a common factor is applied to the real 
and imaginary parts. Turning to > 0, accurate experimental annihilation cross sections 
at low energies are scarce in this mass range, and we only note that the calculated cross 
sections for Ne are in reasonably good agreement with the data 0,^], as can be seen from 
Table |. 

Having gained confidence in the applicability of the density-folded optical potential (F) 
for nuclei, we use it to calculate p reaction cross sections over the periodic table for a range 
of energies, in order to study the systematics of the dependence on Z, on A and on the 
energy. Figure ^ shows, in the upper part, the calculated reaction cross sections for p on 
Ne, Ca, Zr, Sn and Pb at 37.6, 57 and 106.6 MeV/c. It is seen that the range of values is 
very broad; for a given energy the calculated cross sections vary by a factor of between 6 
and 10. For a given nucleus the cross sections vary with energy by a factor between 3.2 and 
5.7. The mechanisms causing these variations are obviously of interest. 

For strongly absorbed particles such as low energy antiprotons, disregarding the Coulomb 
potential, the total reaction cross section may be approximated by ttR"^ where R is the radius 
of the nucleus. The underlying assumption is that antiprotons with impact parameter less 
than R, or equivalently with orbital angular momentum up to /max where semiclassically 
/max + 1/2 = kR, are totally absorbed. However, due to the focussing effect of the attractive 
Coulomb potential, particles with impact parameters larger than R also interact with the 
nucleus, thus causing the cross section to increase. The low energy Coulomb problem is 
well described by the semiclassical approach |]10| as recognized a long time ago by Blair in 
connection with a particle reactions on nuclei |2^. Assuming that total absorption occurs 
in all partial waves for which the distance of closest approach is smaller than R, one gets 
the following relation between the Coulomb-modified /max and R: 

{im..+lr^{kRni+^) . (21) 

At very low energies, 21] >> kR, and therefore /max >> kR due to the focussing effect of V^. 
The total reaction cross section in the strong absorption limit is then given by 
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^ ^ -H 1) ^ i)^ « + H) ^ .fi^d + . (22) 

where and are the laboratory and cm. wave numbers, respectively. The second 
term within the brackets represents the Coulomb focusing effect and at very low energies 
it becomes dominant 0], thus leading to an A^^^Z dependence of the cross section if i? = 
rgA^/^. At high energies the usual strong absorption value ttR^ is obtained. In order to use 
this expression one needs to define an equivalent radius R that will properly represent the 
p-nucleus interaction at the particular energy. We proceed to define such a radius. 
Expanding the p-nucleus scattering wave function in partial waves, 

^{r) = J2{2l + l)^^Pi{cose) , (23) 

we define an average radius < r; > for the I partial wave 

JrW{r)\iJi{r)\^dr 
' JW{r)\^i{r)\^dr ^ ^ 

with W{r) = — Im Vopt(r). Then, bearing in mind the partial- wave expansion of the reaction 
cross section Eq. (p^, we define an average radius as follows: 

E(2/ + i)T, ■ ^^^^ 

Expanding alternatively the reaction cross section Eq. (H) in partial waves: 

2 



aR = —jw{r)A7^Y.i'^l + l)\iPi{r)\'dr , (26) 



T, = ^-^ lW{r)\Ur)?dr , (27) 



one finds that 



which combined with Eq. (^) leads to an equivalent expression for < r >: 

SW{rmv)\^rdv 

!W{r)\ij{v)\Hr ■ ^ 

Note that the average radius turns out to be slightly energy dependent. At pl = 37.6 
MeV/c, the lowest p incident momentum for which pp annihilation has been measured [IH, 
the calculated values of < r > appropriate to the density-folded optical potential (F) can 
be parameterized to better than 1% by 

< r >= 1.840 + 1.120A^/^ fm . (29) 

In fact, this expression holds to better than 2% for the whole range of very low energies up 
to PL ~ 100 MeV/c. 

In the lower part of Fig. |^ are plotted the ratios of the calculated cross sections to 
the semiclassical expression (|2^ ) where for R we took either the radius of a rigid sphere 
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having the same average radius as the optical potential (F), namely i? = | < r >, or 
just R =< r > which could be more appropriate due to the substantial attenuation of the 
wavefunction which is involved in the definition of < r >. It is seen that these ratios vary 
remarkably little, by less than 5% for the lower band, when the cross sections vary by factors 
between 3 and 10. At the lower energy the Coulomb focusing is dominant but at the higher 
energies covered here both terms of Eq. (^) are important. It may therefore be concluded 
that the dependence of the calculated reaction cross sections on the three parameters {Z, A 
and energy) is almost exclusively given by Eq. (^21). 

Another consequence of the strongly absorptive p-nucleus optical potential is the satu- 
ration of the physical observables which are induced by the absorptivity. Thus, for < 0, 
widths of p atomic levels saturate as function of the magnitude of Im V^pt and, for 
E > 0, total reaction cross sections also saturate with it [|]. As argued in our earlier work 
0], the saturation property of for any given nucleus also implies that cjr does not rise 
with Z and A as fast as it would have risen if the p absorption were considerably weaker 
and therefore could be treated perturbatively. This is demonstrated here differently, in Fig. 
^, where the saturation factor S defined as 

g ^ gfl(Imfeo)/Im&o 

aR{10-Hmbo)/lO-^lmbo ^ ' 

is plotted at pi = 57 MeV/c (1.73 MeV) as function of A, using the density-folded optical 
potential (F) for A > 10, (F') for ^He (see below in Sect. |V|) and (G) for H. For Im bo as 
small as 10"'^ of its nominal value (of 1.2 fm for potential (F)), the calculated an is very 
nearly linear in Im bo. If the linear rise of were sustained up to the nominal value, the 
saturation factor S would have assumed the value S = 1. However, for Im bo = 1.2 fm for 
potential (F), the calculation is strongly nonperturbative, as witnessed in the figure by the 
small values of S relative to 1, about 0.1 and less beginning with Ne. The saturation of 
has two aspects to it: (i) that S depends very weakly on Im bo in the nonperturbative regime 
(not shown in the figure); and (ii) that S monotonically decreases with A, as is shown in the 
figure, both for antiprotons (solid curve) and for antineutrons (dashed curve). The values 
of S for antiprotons are lower than for antineutrons due to the Coulomb focussing effect 
which for antiprotons acts in two seemingly opposing ways: (i) enhancing the reaction cross 
section in the negatively charged case compared to the uncharged case, but also (ii) expelling 
more effectively the p wavefunction from the nuclear region, so that becomes more 
strongly suppressed than for antineutrons. This latter effect gets weaker with increasing 
energy; at = 6.04 MeV, Sp is larger by 20 - 25 % than at 1.73 MeV, whereas Sn is 
smaller than at 1.73 MeV. 



V. THE ANTIPROTON - HELIUM SYSTEM 

In this section we study the p - '^'^He system at very low energies, within the optical model 
approach, in order to test the scattering-length approximation (SLA) which was successfully 



applied to the pp system |28| and has recently been extended to the pd and p - ^He systems 
1^. ^He is one of the very few nuclei for which data on the strong interaction of p are 
available for both the very low energy scattering regime and the atomic negative energy 
regime. There are two data points for the total annihilation cross section below 3 MeV in 
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the positive energy regime M , and three data points for the atomic 2p level shift and width 



plus the width of the 3d level |29[ above —20 keV in the negative energy regime. We chose to 
test the ability of the optical model to reproduce these data together with the corresponding 
set of atomic shift and widths data for ^He p9[, although the optical model approach is not 



very well suited for studying such light nuclear targets. 

Since potential (F) does not reproduce very well the p-nuclear observables for A < 10, 
we searched over its parameters to fit the p -^'^He data mentioned above. As was already 
encountered in the fitting procedure described in our earlier work 0] for potential (a) of 
that reference, when trying to fit the atomic 2p and 3c? levels simultaneously the calculated 
width of the 2p state in '^He came out too small. This has been noted by other researchers 



before, as summarized in Ref. |]29|. We therefore dropped out the 3d widths data from the 
fitting procedure, partly on the grounds that the rf-wave contribution to the very low energy p 
annihilation cross sections, which are also being explored here, is almost negligible compared 
to the dominant s and p waves (see below). The resulting density-folded optical potential, 
referred to as (F'), has the following parameters: = 1.8 fm for the range parameter of the 
two-body Gaussian interaction folded in with the He densities, and a strength parameter 
bo = —0.26 -|- 22.07 fm in the notation of Eq. (^. The optical potential (F') is more 
absorptive than potential (F) and, again, its real part plays only a minor role. The 2p 
level shifts and widths, as well as the p total annihilation cross sections calculated using 
the optical potential (F') are shown in Table ||, where the measured values are also given, 
including a recent report of p annihilation on ^He Clearly, potential (F') describes well 



these data. We note that omitting the ^He data from the fit hardly changes the resultant 
(F') which describes very well the "^He data. In particular, it was not possible to decrease, 
within such a restricted fit, the range parameter qq from its relatively large value given 
above. 

Protasov et al. have recently fitted the two low-energy p - ^He total annihilation cross 



sections listed in the table, using the SLA expressions for s, p and d waves ||2H]. The input p- 



and ci-wave Coulomb-modified scattering 'lengths' were derived using the Trueman formula 



1^ from the 2p shift and 2p and 3d widths. The s-wave Coulomb-modified scattering length 



a'o'^^ was left as a fitting parameter, since the p atomic Is level shift and width in He are 
not known experimentally. Assuming a value of 1.0 ± 0.5 fm for Re cto*'^'', these authors were 
able to fit the annihilation cross sections with the following value for the imaginary part of 
the s-wave Coulomb modified scattering length: 

Im 4'^^ = -0.36 ± 0.03(stat.)lS:n(syst.) fm . (31) 

The quality of the fit, as is evident from Fig. 3 in Ref. is similar to ours (see Table |I| 
above), with the calculated lower (higher) energy cross section somewhat above (below) the 
mean value of the measurement. In contrast, our p - ^He potential (F'), which is also fitted 
to essentially the same data set, upon using Eq. (|T^) yields the following values: 

Re ai'^^ = 1.851 fm , Im af,'"^ = -0.630 fm . (32) 

It is clear from this example that the 'model independent' determination claimed by Protasov 
et al. for a^^'^^ is violated by our specific example, suggesting that, contrary to their claim, 
it is not a model independent determination. In order to study the origin of the above 
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discrepancy, we show in Table |TT| the partial wave contributions to the calculated cross 
section at pi = 57 MeV/c. The most important contributions are due to the s- and p- 
waves, for which the two calculations give quite different results in spite of agreeing well 
with respect to their sum. In particular, the p-wave contributions are quite different from 
each other, although sharing practically the same value of the Coulomb-modified p-wave 
'scattering length' (more traditionally called 'scattering volume') a^i^\ Whereas at —20 
keV, the I = 1 p - ^He dynamics is well determined by a^*'^^ alone, over the energy range of 1 
- 3 MeV corresponding to the annihilation measurements, it depends on more than just this 
'scattering length'. The effective range term, and perhaps higher order terms in the effective 
range expansion at low energies, become equally important. Indeed, we have verified for 
potential (F') that the variation of the Coulomb-modified / = 1 scattering phase shift is not 
reproduced in this energy range by specifying the 'scattering length' alone. Once the p-wave 
contributions to the total annihilation cross section differ by as much as is observed in the 
table, the essentially fitted s-wave contribution in the calculation of Ref. must also differ, 
and hence our prediction for the Coulomb- modified s-wave scattering length is necessarily 
different from that of Ref. 0. Such a difference would not occur for the pp system, which 



at the appropriate low energies is largely controlled by s waves ||28[ . 

Protasov et al. also attached particular significance to the observation that, for the 
Coulomb-modified s-wave scattering length, "the absolute value of the scattering length 
seems to be a decreasing function of the atomic weight" since "a naive geometrical picture 
of p-nucleus annihilation would suggest a value of the p-nucleus scattering length increasing 
with the nuclear size." Indeed, their expectation holds for the plain strong-interaction (not 
the Coulomb-modified) scattering lengths, as borne out by the calculations of Ref. which 
are updated in the next Section [V^. However, the Coulomb-modified scattering lengths do 
not show such a clear geometrical picture, as is also discussed in Section 0. 



VI. ANTIPROTON - NUCLEUS SCATTERING LENGTHS 



In an earlier paper |^ , hadron nucleus scattering lengths were derived from exotic atom 
data for kaons, antiprotons and sigma hyperons. For the purposes of the present paper, 
we have updated the previous calculations for antiprotons, using more recent experimental 
measurements of strong interaction shifts and widths. 

For exotic atoms with Z > 1, strong interaction measurements are only available for 
angular momentum states with / > 0. As a result it is necessary to use an optical model 
approach as an intermediate step in deriving the s-wave Coulomb-modified scattering length 
or the purely strong-interaction scattering length Oq from the measured strong inter- 
action effects for I > atomic states. In this method, the optical potential given by Eq. 
(0) was used to fit the strong interaction shift and width values for individual nuclei by 
adjusting the real and imaginary parts of the optical strength parameter bo. The s-wave 
scattering length oq was then calculated by solving the Klein-Gordon equation for I = 
with the optical potential Kpt, but without the Coulomb potential, at an energy (1 keV) 
close to threshold. Further details of the method are given in the earlier paper . 

The experimental shift and width measurements were taken from the published literature 
and cover the available target elements from C to Pr inclusive, omitting Yb since the nucleus 



is deformed. The data set used is discussed and tabulated in []8|,|32[|. For the present work. 
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single particle distributions were used P,|32| for p(r) as these are expected to be more 



appropriate for the analysis of antiproton data. In the previous analysis macroscopic 
density distributions were used. 

The results of the present analysis are shown in Fig. ^ and are qualitatively similar to 
those obtained earlier [^], where an opposite sign convention was used for ao (there denoted 
as). Again the real part of the scattering length shows an approximate A^^^ dependence, and 
the imaginary part is constant. This is to be expected on the basis of a simple model based 
on a strongly absorbing square well potential A least squares fit to these scattering 

length values gives: 

Re ao = (1.54 ± 0.03)^°-^"^° °°^ fm , Im ao = -1.00 ± 0.04 fm . (33) 



These best fit parameter values are in excellent agreement with those obtained earlier ||31|| . 
We note that the magnitude of Im oq is considerably larger than expected for a sharp-edge 
potential, resulting mainly from the diffuseness of the potential ||3T| , |33[| . 

For completeness, values for ao for nuclei with A < 10 are also shown in Fig. ^. The 
values for ^Li and ^Li were determined using the method described above from the strong 
interaction shift and width measurements of Poth et al. The values for ^He and ^He 
were calculated using potential (F') described in Sect. |V|. The value for hydrogen was 



obtained with potential (G) described in Sect. |T|. The real parts of ao for these light nuclei 



are seen to deviate from the predictions of the best fit parameters for A > 10 described 
earlier, whilst the imaginary parts remain relatively independent of A. 

To confirm this behaviour for nuclei with A < 10, the analysis of the strong interaction 
data for C to Pr p|,p2|, together with that for ®Li and ^Li |Q, was repeated fitting all the 
data simultaneously with a single value for the complex parameter Bq. A good fit to the 
strong interaction data is obtained with a of 52 for 46 data points. The s-wave scattering 
lengths obtained in this way are typically a factor of 5 more precise than those in Fig. ^, 
which were obtained using fits to individual nuclei, but they are model dependent since it 
is assumed that the data for all nuclei can be fitted with a single value of bo. However this 
seems to be a reasonable assumption in view of the very good fit to the data. The results 
confirm the deviation in the value of Re ao for ^Li and ^Li from the simple parameterization 
for A > 10 nuclei shown in Fig. ^(top). The absolute values of Im ao of 1.27 ± 0.04 and 
1.16 ± 0.04 fm for ^Li and ''Li respectively, are larger than the mean absolute value of Im 
ao = 1.04 ± 0.01 fm obtained for A > 10. 

To conclude this section, we wish to study the A dependence of the p-nucleus Coulomb- 
modified s-wave scattering length ao**'^'' for a given potential, here chosen as potential (F). 

(sc) 

The calculated values of ao for light and medium-weight nuclei are given in the first two 
rows of Table The behaviour of these scattering lengths with A does not follow a clear 
geometrical picture. The real part at first increases, it then decreases very quickly and flips 
sign, becoming negative already for ^^C. Another sign flip occurs somewhere between Ne 
and Ca. This oscillatory behaviour of Re ag'^^'* is very different from the smooth rise of Re ag 
with A (Eq. ( p^ ) and Fig. The magnitude of the imaginary part at first increases with A 
(rather than decreasing according to the speculation of Ref. 0), reaching a maximum value 

(sc) 

somewhere near A ~ 11 — 12 approximately where Re ao vanishes, it then decreases to a 
minimum value somewhere between Ne and Ca. This behaviour of Im ao'^'^'' is very different 
from the approximate constancy established above for Im ao- 
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The behaviour of function of A can be quahtatively understood using the approx- 

imate relation Eq. ( pOD between oq and a'o'^^ which assumes that R <^ as, where R is the 
range of the short-ranged potential giving rise to a^, and is the Bohr radius. Neglecting 
(Im a^y relative to (Re Oq)^, approximating R by Re Uq ~ 1.5A^^^ fm and keeping only the 
logarithmic term within the brackets in Eq. (|20|), the condition for Re a'^'^^ = is given 
by Int ^ —t, where t = a^'*/(3A^/^) with a^'* = 57.64 fm for the Bohr radius of pp and 
assuming Z = A/ 2. This condition is satisfied for A ^ lA for which R = 3.62 fm according 
to the above approximations, so that the condition of applicability R <^ as = 4.41 fm is 
still (although barely) satisfied. Relaxing the numerically unnecessary approximations of 
neglecting Im oq and (1 — 27), the vanishing of Re Oq'^'' occurs for A ^ 13.5, where Im Qq'^^ 
attains a substantial maximum. Thus, Eq. ( pOD explains well, within its limited range of 
validity, the main features of Re a'^'^^ and Im al^'^^ from Table The 'non-geometric' 
dependence of a^^^ on A is due to the logarithmic term in Eq. (pOf ) which arises from the 



unavoidable logarithmic behaviour of the Coulomb wavefunctions ||1CI|| . 

The s-wave Coulomb-modified scattering lengths Qq'^'' correspond to the short-ranged 
potential which, strictly speaking, is superposed on a point Coulomb potential yj^"*"*-* = 
—Za/r. This short-ranged potential consists then, in addition to the strong-interaction 
optical potential Vopt, of the short-ranged Coulomb term V^(r) — \/^p°*"*) (r) . In practice we 
have also included a vacuum-polarization potential in the latter short-ranged contribution. 
The Coulomb-modified phase shifts Si^'^^ due to the short-ranged part of the interaction, 
and the related scattering 'lengths' a^^^^ considered in Sect. ||, therefore include also effects 
which are not entirely due to the strong interactions. To be more precise, the superscript 
(sc) should be replaced by {st + fs;pc), that is strong + finite size with respect to point 
Coulomb. It would be more fitting, perhaps, to define the Coulomb-modified phase shifts 
^(stjs+pc)^ due exclusively to the strong interactions with respect to the overall Coulomb 
potential V^. Similarly, one could also define the 'finite-size' Coulomb-modified phase shifts 
^Y^^ short-ranged non strong- interaction potential with respect to V'J^"*"*-'. It 
is straightforward to show that 

^(sc) _ ^(st+fsypc) _ ^{st;fs+pc) _|_ ^(/s;pc) ^g^-^ 

Using Eq. (^) and similarly tailored low-energy effective-range expansions, Eq. (0) leads 
to the following relation between the corresponding scattering lengths: 



(sc) (fs;pc) 

~{sc) ^ (st-Js+pc) _ ap - Q-0 /oc\ 

1 + {2Tc/aB) % 

where a^^'^^ = a'^^'^'^^'^^^ is the same Coulomb-modified short-ranged scattering length hith- 
erto considered and tabulated in the first two rows of Table a^/^'^^^ is due to the finite 
size Coulomb (plus vacuum polarization) effect, and Oq'^^ = a\^^''^^^'^^'' on the l.h.s. is the 
intrinsically strong-interaction Coulomb modified scattering length. We have also tabulated 
these newly defined scattering lengths in Table The real finite size scattering length 
j^jggg modestly with Z, becoming significant for ^^C where it leads to a change from 
—0.6 fm for Re oIq^'' to 0.9 fm for Re ciq'^\ The induced change for Im Oq'*^'' into Im Oq^^ al- 
ready becomes significant for ^Be. However, the general trend of Oq^'^ with A is qualitatively 
similar to that described and discussed earlier for cto*"^'*. 
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VII. SUMMARY AND CONCLUSIONS 



We have shown in this work that a unified optical-model description of low-energy an- 
tiproton interactions can be successfully achieved using a density-folded p optical potential 
Vopt ^ V * p. As shown in Sect. |Tl|, the potential v gives an excellent fit to the measured 
pp annihilation cross sections at low energies {pl < 200 MeV/c) and to the Is and 2p 
spin-averaged pB. level shifts and widths, but is far from being unique; other combinations 
of range and of strength parameters have been derived for the pp system for low as well 
as for higher energies ||I9|J20[| , all of which yield a highly absorptive potential. It would be 



useful to measure also pp elastic scattering differential cross sections at these low energies, 
including the Coulomb-nuclear interference angular region, in order to attempt resolving the 
ambiguity in the determination of v, particularly its less well determined real part. 

The optical potential Vopt reproduces satisfactorily the p atomic level shifts and widths 
across the periodic table for A > 10 as well as the few annihilation cross sections measured 
on Ne. This simple folding model does not work well for the very light species of He and 
Li, possibly due to the spin and isospin averaging and to other unjustified approximations 
inherent in the construction of optical potentials for these relatively small values of A. 
However, for a limited range of small- A values, as shown in Sect. for the He isotopes, 
an energy-independent optical potential can be fitted to both the measured annihilation 
cross sections a.t E > and the atomic shift and width data for E < 0. Our detailed 
analysis of the p—'^B.e system was primarily intended to demonstrate the inapplicability of 
treating it by the scattering length approximation 0. As a byproduct, we have discussed in 
Sect. the systematics of several sets of p - nucleus s-wave scattering lengths, the purely 
strong-interaction set for which we confirmed and updated the discussion of Ref. [^|, and 
the Coulomb- modified set which defies any geometrical interpretation. 

Fairly accurate reproduction of the optical potential calculations of was achieved 
in Sect. |I^ by extending the black-disk strong- absorption model to account for the 
Coulomb focussing effect of negatively charged projectiles at very low energies {pl < 100 
MeV/c). As the number of additional partial waves contributing to the p reaction cross sec- 
tion due to the Coulomb focusing effect increases substantially upon decreasing the incoming 
energy, the semiclassical description becomes valid, resulting in the simple closed-form 'ge- 
ometrical' expression (^) which embodies the overall dependence of on Z, A and on 
E, provided the black-disk radius is properly handled (see Eq. (pQ])). Figure ^(bottom) 
demonstrates this success for very low p energies over the periodic table. 

A final comment is due to the apparent success of using a density folded optical potential 
in terms of a potential v that describes well the pp system. Deloff and Law advocated and 
used a. V * p folding model for kaonic atoms [§^,0 and for p atoms [^. This model was 
subsequently tested by Batty for kaonic atoms and for antiprotonic and sigma atoms 
and moderately good agreement with the data was found. However, the underlying 



potentials v were not fitted to a comprehensive set oi Z = I data. The present work 
provides a step forward for antiprotons, in that a complete body of pp data was fitted to 
determine v. This has to be supplemented in due course by pn data, particularly in order 
to understand the renormalization factor of 2/3 encountered in obtaining the density- folded 
optical potential (F) from the v potential (G). A multiple-scattering justification of the v* p 
folding model was given by Green and Wycech |^ for antiprotons as due to the strong 
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absorptivity, and was confirmed in detailed calculations by other authors pl] , ^ . 

Clearly, more data on p nuclear interactions at low energies are needed to develop quan- 
titatively further the unified optical-model approach put forward in the present work. Such 
data should consist of elastic and charge exchange scattering, as well as annihilation cross 
sections on nuclei. In particular, data on the very light nuclei, for A < 10, would provide 
valuable information on how to bridge the gap, prevailing at present within the density- 
folded optical model approach, between hydrogen and the A > 10 nuclei concerning the 
applicability of the same (smoothly dependent on A) p optical potential. 

CJB wishes to thank the Hebrew University for support for a visit during which this 
work was started. 

This research was partially supported by the Israel Science Foundation. 
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FIG. 1. Top: experimental annihilation cross sections [1,16] for the pp system and calcu- 
lated values using the (G) and (SW) potentials. Bottom: ratios of the real part of the pp 
Coulomb-modified forward scattering amplitude to its imaginary part for the same potentials. 
Real to imaginary ratios for the zero-energy scattering length are also shown. 
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FIG. 2. Top: calculated p-nucleus reaction cross sections for three incoming momenta using 
the density- folded potential (F). Bottom: ratios of the above cross sections to the semiclassical 



expression Eq. (pi). Upper band for R =< r >, lower band for i? = | < r >. 
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FIG. 3. The saturation factor S of Eq. (|30|) at pi = 57 MeV/c, calculated using the den- 
sity-folded optical potential (F) for A > 10, (F') for A = 4 and (G) for ^ = 1, as function of A for 
antiprotons (solid line) and for antineutrons (dashed line). 
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FIG. 4. Real (top) and Imaginary (bottom) parts of the s-wave p scattering length oq calculated 
by fitting to p atomic data. The straight lines are a best fit to the calculated values for A > 10, 
see Eq. (H). 
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TABLES 



TABLE L Measured and calculated annihilation cross sections for p on Ne. Calculations were 
made with potential (F). 

PL (MeV/c) dimf^ (mb) 

57 2210±1105" 2524 

192.8 956±47^ 1016 



(Tann (mb) 



Ref. [i 
^ ReL [|6| 



TABLE IL Measured and calculated observables for the p— ^'^He system. Shifts (e) and widths 
(r) are in eV, cross sections are in mb. Calculations were made with potential (F'). 







e2p 


^2p 


aann (47.0 MeV/c) 


aann (55.0 MeV/c) 


aann (70.4 MeV/c) 


^He 


calc. 


-12 


33 




1038 




^He 


exp. 


-17 ±4°^ 


25 ±9°^ 




1850±700^ 




^He 


calc. 


-19 


42 


1116 




810 


^He 


exp. 


-18 ±2°^ 


45 ±5°^ 


979 ± 145'' 




827 ± 38^" 



Ref. [H] 
^ Ref. 
^ Ref. M 



TABLE III. Partial wave contributions (in mb) to the calculated p - "^He total annihilation 
cross section at pi = 57 MeV/c. 





/ = 1 = 1 


/ = 2 


/ = 3 sum 


Protasov et al.*^ 


280.3 652.5 


16.2 


949 


potential (F') 


395.8 500.3 


49.8 


1.5 949 


experiment'' 






915 ±32 


« Ref. 








" ReL § 








TABLE IV. Coulomb-modified s-wave p scatttering 


lengths (in 


fm) calculated for potential (F). 


target ^He 


^Li lOg 


12c 


1^0 Ne Ca 



Re a, 
Im a, 





(sc) 




(fs;pc) 

Re a, 
Im a, 



(sc) 


(sc) 



1.562 

-0.590 
0.067 
1.479 

-0.576 



1.967 
-1.012 
0.148 
1.703 
-0.853 



2.169 
-1.792 
0.170 
1.809 
-1.183 



1.401 
-2.853 
0.183 
1.633 
-1.612 



-0.626 
-2.631 
0.176 
0.948 
-2.373 



-0.645 
-0.551 
0.252 
-0.654 
-1.414 



-0.216 
-0.241 
0.307 
-0.485 
-0.556 



0.134 
-0.391 
0.835 
0.004 
-0.121 
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